Controlled absorption and all-optical diode action due to collisions of self-induced 

transparency solitons 



(N 

o 

U 

Oh' 
< 
O 



o 

H— > . 

Or 

q 

o ■ 

• i-H , 

Oh' 



> 

o 

(N 

O 

(N 



X 



Denis V. Novitsk^ 

B. /. Stepanov Institute of Physics, National Academy of Sciences of Belarus, 
Nezavisimosti Avenue 68, BY-220072 Minsk, Belarus 
(Dated: March 2, 2013) 

We study inelastic collisions of counter-propagating self-induced transparency solitons in a homo- 
geneously broadened two-level medium. The energy of the pulse can be almost totally absorbed in 
the medium due to asymmetric collision with a properly chosen control pulse. The medium state 
thus prepared demonstrates the property of an all-optical diode which transmits pulses from one di- 
rection and blocks from another. The saturation process of a controlled absorption effect, local-field 
correction influence, and the parameter ranges for the diode action are studied as well. 

PACS numbers: 42.50.Md, 42.65.Sf, 42.65.Tg 



I. INTRODUCTION 

The problem of creation of all-optical diodes has been 
actively studied during recent years. Such systems, which 
transmit electromagnetic radiation only in one direction 
and block it in another, are usually based on asymmetri- 
cally constructed waveguides [l| or nonlinear photonic 
structures (see, for example, 0t3)- The use of left- 
handed materials @ and asymmetrically absorbing sys- 
tems Q have been studied, among other possibilities. In 
other words, the all-optical diode action is commonly a 
consequence of the proper ordering of optical elements 
on microscopic scale (of the order of a light wavelength) . 

In this paper we propose a fundamentally different 
scheme of all-optical diode action. It is entirely based 
on nonlinear dynamics of the self-induced transparency 
(SIT) solitons in a homogeneously broadened, dense two- 
level medium. SIT is a well-known effect [8, 9] , resulting 
in formation of stationary pulses (solitons) which propa- 
gate without change in their form. Our approach deals 
with collisions of counter-propagating SIT solitons. The 
effects of such collisions were partially studied in our pre- 
vious work jlc} . It was shown that under certain condi- 
tions, one of the colliding solitons can be almost entirely 
absorbed by the medium. We call this phenomenon the 
controlled absorption of the SIT soliton. 

A large density of the two-level medium allows one to 
obtain the resulting effects on shorter distances. It is 
also useful for reducing the scope of the numerical cal- 
culations. The local-field effects need to be taken into 
account in this case. In Ref. 13] it was shown that 
they do not significantly affect the dynamics of a pulse 
with duration of femtosecond order. However, the local- 
field correction is taken into account here for generality; 
moreover, it is possible that it can somehow influence the 
dynamics on relatively long distances, as in our consid- 
eration. Therefore, we discuss this problem further as 
well. 



The paper consists of the following sections. In Sec. 
HT1 we consider the effect of controlled absorption of the 
SIT soliton in some detail. Section HTT1 is devoted to an 
explanation of our mechanism of all-optical diode action. 
In the Sec. [IV] we discuss the influence of pulse launch 
time, relaxation time, and medium density (as well as 
local field) on the demonstration of our effects. Finally, 
in Sec. |V]we make some concluding remarks. 



II. CONTROLLED ABSORPTION 

Propagation of light pulse in the dense two-level 
medium is governed by the system of semiclassical 
Maxwell-Bloch equations for population difference W, 
microscopic polarization R, and dimensionless Rabi fre- 
quency n' — Q/uj = (iJ,/hu))E (i.e., the normalized elec- 
tric field amplitude) [ll|.[l^: 
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where t = cut and £ = kz are dimensionless time and dis- 
tance, respectively; S = Auj/uj the normalized detuning 
of the field carrier (central) frequency u) from the atomic 
resonance; 71 = (wTjJ^ 1 and 72 = (o^) -1 the rates 
of longitudinal and transverse relaxation, respectively; 
k = lo/c the wavenumber; c the light speed in vacuum; 
and e = lol/uj = 4ir fj 2 C / 3hui is the normalized Lorentz 
frequency, where C is the concentration of the two-level 
atoms, and /1 the component of transition dipole mo- 
ment parallel to the polarization vector of the electric 
field. Here we assume without loss of generality that 
the background dielectric permittivity of the medium is 
unity (two-level atoms in vacuum) . Note that in Eq. ([3]) 
we do not use the slowly varying envelope approximation 
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(SVEA), while in Eq. ([T]) the term eW is responsible for 
the so-called local field correction (near dipole-dipole in- 
teraction) in a dense medium. The numerical approach 
used in this paper to solve Eqs. (d])-© can ^ e f° un d in 
Ref. QJ]. 

We consider propagation of ultrashort pulses with 
Gaussian shape ft = ft p exp(—t 2 /2tp), where t p is the 
pulse duration. The amplitude of the pulses is mea- 
sured in the units of the characteristic Rabi frequency 
ilo = ^/2n/2t p , which corresponds to the so-called 2ir 
pulse. In our calculations the values T\ = 1 ns and 
T2 = 0.1 ns are taken so that femtosecond pulses are 
in the regime of coherent interaction with the resonant 
medium. We use the following parameters of calculations 
which hold true throughout the paper: lol = 10 11 s , 
6 = (exact resonance), A = 2irc/ui = 0.5 /im; the thick- 
ness of the medium L = 1000A. All the pulses have the 
same duration t p = 50 fs, which means that the only pa- 
rameter governing pulse dynamics is its amplitude. One 
can estimate that a 2tt pulse with such duration should 
have peak intensity of about 10 MW/cm 2 in the case of 
atomic dipole moments /j~ ID. The required concentra- 
tion in this case is C ~ 10 19 cm -3 . However, the density 
and the peak intensity can be significantly reduced if the 
media with higher dipole moments is used (for example, 
a collection of quantum dots). 

We start with the scheme of two ultrashort (coherent) 
counter-propagating pulses: the forward-propagating 
(FP) pulse with the initial amplitude f2 p i and the 
backward-propagating (BP) one with f2 P 2 (see Fig. [T|). If 
their amplitudes are large enough, both pulses form self- 
induced transparency (SIT) solitons as they propagate 
inside the medium. Such two counter-propagating soli- 
tons meet somewhere inside the medium and interact. It 
is known that, in contrast to the pair of co-propagating 
solitons, this interaction is inelastic Mathematically, 
this means that two counter-propagating pulses do not 
form a stationary solution of the equations of motion. 
From the physical point of view, this fact can be ex- 
plained with the simple argumentation as follows. For 
the co-propagating solitons, we always have a situation 
when the light energy absorbed at the leading edge of 
the pulse is released at its trailing edge. This is not the 
case for the counter-propagating solitons: at the point 
of collision, the medium excited at the leading edge of 
the first pulse interacts with the second pulse even be- 
fore the trailing edge arrives. As a result, the collision of 
the pulses leads to the partial absorption of their energy. 

It was recently shown [l(| that, if = fio (27r-pulse) 
and P 2 = 1.5f2o (37r-pulse), the FP soliton gets entirely 
absorbed, while the BP one appears unperturbed at the 
output. In other words, the first pulse is blocked by 
the second one. For simplicity, we call this situation the 
controlled absorption of light energy inside the two-level 
medium. This energy leaves the medium in a long run as 
a result of fluorescence. We should also emphasize that 
this effect occurs only for the asymmetric collision, i.e. 
when the amplitudes of both pulses are not the same. 
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FIG. 1. The scheme of collision of two counter-propagating 
(FP and BP) pulses in a two-level medium. 
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FIG. 2. (Color online) The resulting intensities after the 
collision: (a) collision of two pulses with fi p i = Qo and 
Q, P 2 = 1.5fio; (b) the third (probe) FP pulse Q P 3 = fio is 
added; (c) the probe FP pulse has Q p z = 1.5f2o; (d) the probe 
pulse with f2 P 3 = flo is BP. The first and second pulses start 
at the instant t = 0, while the third one starts at t — 100i„. 



The resulting profiles of the FP and BP pulses with ini- 
tial amplitudes = Qq and Q P 2 = 1.5f2o are demon- 
strated in Fig. [2a). These profiles are recorded at the 
opposite edges of the medium: according to Fig. [TJ it 
is right edge for the FP pulse and the left edge for the 
BP one. One can easily see that the FP pulse is absent 
at the output (its energy is mostly absorbed), while the 
BP pulse appears at the other end of the medium. Note 
that this pulse (we can call it the control pulse) can be 
described by the usual 2-7r-soliton envelope and is signif- 
icantly compressed in comparison with the initial pulse 
due to the transient process [l(|. There are also low- 
intensity tails connected with radiation of absorbed light 
at both outputs of the medium. 
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III. DIODE ACTION 

Now one can ask a question: what happens when the 
third (or probe) pulse with the amplitude tt P 3 enters the 
medium? There are three possible answers: (i) it will 
be trapped, (ii) it will pass through the medium, and 
(iii) it will not only pass but also retrieve the previously 
absorbed pulse from the medium. Further we try to in- 
vestigate this question for different values of S7 P 3 . If the 
possibility (iii) can be realized, one can talk about optical 
memory 

Let us add a probe pulse with amplitude f2 P 3 which en- 
ters the medium after the interval t = 100t p when there 
is no propagating solitons inside the medium. There are 
two alternatives: the probe pulse is FP or BP. The case 
of FP probe pulse with f2 P 3 = fig is shown in Fig. [2jb). 
It is seen that interaction of this pulse with the medium 
(which stores the energy of the previous FP pulse) pre- 
vents its appearance at the output. Only tails have some- 
what increased intensity. As opposed to this case, the 
pulse tt p 3 = 1 .5ilo is happily transmitted through the 
medium though with slightly lower peak intensity than 
the identical BP pulse [Fig. 0c)]. This means that trans- 
mittance of the soliton depends on its intensity. More- 
over, we have an interesting deviation from the area the- 
orem: the pulse with initial area 2ir gets absorbed, while 
the one with 37r forms a 27r-soliton. 

Now, if the probe pulse is BP and fi P 3 = f2g, then it 
appears at the output [Fig. 0cl)], in sharp contrast with 
the FP case [Fig. HJb)]. In other words, the collision 
of two counter-propagating pulses switches the medium 
into such a state that it transmits the probe pulse inci- 
dent from one side and blocks from the other. We can call 
these effect the all-optical diode action. The reason for 
this effect can be understood if we recall that the colliding 
pulses are not identical, i.e. the collision is asymmetric. 
As a consequence, this asymmetric collision "prepares" 
the asymmetric state of the medium, i.e. the diode ac- 
tion entirely depends on light-matter interaction, but not 
on the properly chosen geometry of the problem. There- 
fore, this scheme is fundamentally different from usual 
proposals of all-optical diodes based on asymmetrically 
constructed systems. 

To study this asymmetry in more detail and to clarify 
the conditions for the diode action, we investigate the 
dependence of absorbed energy on the amplitude of the 
probe pulse f2 P 3 (see Fig. [3]). Absorbed energy is cal- 
culated as a part of total electromagnetic energy which 
remains inside the medium after some period (namely, 
t = 200t p ) from the start of the first two pulses. The 
curves for both (FP and BP) cases are situated between 
two asymptotes. The first one characterizes the fraction 
of the total energy of all three pulses which accounts for 
the first FP pulse. This asymptote is defined as the func- 
tion ft(£lpa) 
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p i = Hq, f2 P 2 = 1.5f2o an d 
£! p 3 = Qp^/flQ. It describes the absorbed energy in 
the ideal case when only the energy of the first pulse is 
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FIG. 3. (Color online) The dependence of absorbed energy 
(as a fraction of the total energy) on the amplitude S1 P 3 of the 
probe pulse. The energy is integrated over the time period 
t = 200t B . 
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FIG. 4. (Color online) The dependence of peak intensity and 
propagation time (in units of t p ) on the amplitude Q p $ of the 
probe pulse. For comparison, the results for a single pulse 
with the amplitude Q P 3 are shown. 



trapped in the medium. Similarly, we introduce the sec- 



ond asymptote / 2 (fi p3 ) = (^1+^3)/ (^1+^2+^3) = 



L p3)/ V"pl^"p2^"p3^ 

(1 + Qp 3 )/(3.25 + ^p 3 ), which describes the part of the 
energy carried by the first and third pulses together. It 
corresponds to another ideal situation when both the first 
and the probe pulses are absorbed. 

It is seen in Fig. [3] that at low values of f2 P 3, the ab- 
sorbed energy curves for the FP and BP probe pulses 
virtually coincide and are situated closer to the /2- 
asymptote, i.e. a large part of energy of the third pulse is 
absorbed inside the medium. However, at f2 P 3 > 0.8f2g, 
the curves sharply diverge: the FP probe pulse is still 
mainly trapped, while the BP probe pulse is almost to- 
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FIG. 5. (Color online) The dependence of absorbed, trans- 
mitted (FP), and reflected (BP) energy on the number of 
additional probe pulses with the amplitude Qo- The energy 
is integrated over the time period t = 1000t p . 



tally transmitted (the blue curve with triangles gets very 
close to the /i-asymptote in Fig. [3}. The energy of 
the first pulse is absorbed by the medium in both cases. 
Thus, the most distinction between the FP and BP cases 
is observed exactly in the region f2 P 3 w Hq. This is the 
asymmetry seen in Figs. GJb) and[2jd). At larger values 
of Slp3 , the curves converge again, so that the medium be- 
comes transparent for the FP pulse as well [cf. Fig. HJc)]. 
Finally, for the FP pulses with 51 P 3 ss 2fio, the fraction 
of trapped energy increases again. This allows us to see a 
certain periodicity of absorption, though we will not con- 
sider the case of large intensities (£l P 3 > 2Qq) here. One 
can expect more sophisticated dynamics in that case due 
to the process of pulse splitting into several 2n solitons. 

The distinction between the FP and BP probe pulses 
is easily seen in Fig. |4l which shows the peak intensity 
of the probe pulse at the output versus its amplitude at 
the input. The BP probe pulse behavior almost coincides 
with that of the single pulse (i.e., with O p i = fi P 2 = 0). 
On the contrary the FP probe pulse has significantly 
lower value of peak intensity, which is equal to zero at 
f2 p 3 < 1.5f2o (no soliton at the exit). Note that at high 
input amplitudes, the secondary pulse appears as a re- 
sult of pulse splitting of both the single and BP pulses. 
However, this is not the case for the FP probe pulse, i.e., 
the controlled absorption blocks the splitting of a high- 
intensity pulse in this case. 

The next question to be discussed is as follows: How 
many FP pulses can be absorbed inside the medium? 
Obviously, the number must be limited. To find the an- 
swer, we launch a number of identical FP pulses (with 
f2 p 3 = Qq) in 100ip one after another and examine the 
fraction of absorbed energy as a function of probe pulse 
number n. The asymptote for ideal trapping of all the 
pulses is fs(n) 
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FIG. 6. (Color online) (a) Transmitted radiation profile for 
the nine incident probe pulses with the amplitude Qo- (b) 
Distributions of population difference W inside the medium at 
the instant t — 1000t p for different numbers of probe pulses. 
Population difference W = 1 corresponds to atoms in the 
ground state. 



(1 + 7i)/(3.25 + n). The results of calculations for dif- 
ferent n are presented in Fig. [5j It is seen that the dif- 
ference between the asymptote and the absorption curve 
grows rapidly with n. Moreover, the absorbed energy 
stops increasing already at n > 3. This means that the 
saturation of light trapping takes place, while simultane- 
ously the transmitted energy slowly grows. The increase 
of transmitted energy becomes significant at n > 6, while 
the absorbed part of the energy begins to decrease. One 
can expect that there should be not only the low-intensity 
tails at the output, but also the solitons as well. This ex- 
pectation is proved to be valid by Fig. [H^a): there is no 
transmitted pulse at n = 6, but, starting from n = 7, 
every extra pulse appears at the output. 

Figure [B^b) allows us to trace this saturation process. 
It is seen that the energy of the first absorbed pulse 
(n = 0) is localized inside the medium: medium exci- 
tation near 200 /jm is due to the collision with the BP 
pulse, while the postcollisional pulse is trapped at larger 
distances (mainly, at L < 400 /mi). As we launch addi- 
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FIG. 7. (Color online) The dependence of absorbed energy 
(as a fraction of the total energy) on ts, the launch time of 
the probe FP pulse with amplitude S1 P 3 = Oq- The energy is 
integrated over the time period t — 1000t p . Time is measured 
in units of t D . 



tional (probe) FP pulses, their energy is absorbed in the 
unexcited regions of the medium (n = 3) until the distri- 
bution of excitation becomes practically uniform (n — 6). 
Note that the region L > 400 /im remains almost unex- 
cited, which was predetermined by the first pulse. This 
uniform distribution of population difference corresponds 
to the saturation and practically does not change when 
we add more pulses {n — 9). As a result, the interaction 
of the additional pulses gets weaker for n > 7, so that 
they can pass the medium. These smeared distributions 
also allow an understanding of why the absorbed pulse 
cannot be recovered. Therefore, we cannot talk about 
the storage of light and use the term " controlled absorp- 
tion" . The absorbed energy is radiated by the medium 
dipoles during the interval of the order of the relaxation 
time. 



IV. PARAMETERS REGIONS FOR THE DIODE 
ACTION EXISTENCE 

In the above consideration we adopted the time inter- 
val 100ip between the launch instants of the first two and 
the third pulses. Here we must point out the significance 
of the time interval between the pulses, which must be 
much less than the relaxation times to preserve the coher- 
ence of light-matter interaction. We study this temporal 
effect in Fig. which shows the absorbed part of the 
energy as a function of the launch time £3 of the probe 
(second FP) pulse. The first two pulses were launched 
at the instant t\ = t% = 0. It is seen that at short 
<3 < 200t p , absorbed energy is somewhere between the 
energy of the single FP r2 pulse (0.235) and the energy 
of both FP fin pulses (0.47). This means that at the in- 
stant t = lOOOtp (the final time of energy integration), 
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FIG. 8. (Color online) The dependence of (a) the peak inten- 
sity and (b) propagation time (in units of t p ) on tz, the launch 
time of the probe FP pulse with amplitude Q P 3 — fio- The 
cases of absence of local-field correction (LFC) and relaxation 
are shown as well. The "standard curve" is calculated taking 
into account both LFC and relaxation. 



the medium contains the energy of the first FP pulse 
and partly of the probe pulse. This absorbed energy 
slowly decreases as we increase the launch time t%, and 
at £3 > 300i p it drops below the level of the single FP 
pulse energy. One can suggest the main role of relaxation 
in this behavior. Indeed, the time, say, t — 500t p = 25 
ps, is comparable with the relaxation time T 2 = 100 ps 
adopted in our calculations. To prove this suggestion, 
we performed simulations for the relaxation-free medium 
(Ti = T 2 = 0). The results are presented in Fig. (the 
dashed line with the triangles). It is seen that the curve 
in this case keeps above the level 0.235, i.e., the relax- 
ation does not affect the process of controlled absorption 
of the first FP pulse. However, the energy of the probe 
pulse is not absorbed as well. Nevertheless, as Fig. [5] 
demonstrates, the peak intensity of the pulse transmit- 
ted is strongly suppressed for the relaxation-free medium 
in comparison with the case of the presence of relaxation 
(this last case is shown by the "standard curve" in the 
figure). Thus, relaxation helps to form an intense soliton 
during propagation through the medium prepared by the 
collision of two counter-propagating pulses. 

Figure [8] also clarifies the influence of local-field cor- 
rection [the nonlinear addition to the detuning in Eq. 
(UJ) ] on the effects discussed. It is seen that this in- 
fluence can be considered as negligible for our param- 
eters. This is in full accordance with the conclusion of 
Ref. [13| . Indeed, for the parameters of our calculations 
(ujl = 10 11 s _1 , t p = 50 fs), the ratio of the pulse ampli- 
tude (peak Rabi frequency) to the Lorentz frequency is 
Qq/wl = V^tt ■ 10 2 >> 1, while, for the local- field effects 
to be observed, one needs Qo/^l ~ 1- 

Now let us consider the effect of pulse duration on the 
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duration of the incident pulses. 

Obviously, the dynamics of light depend not only on 
the parameters of the pulses themselves, but also on the 
parameters of the medium, such as its density. In other 
words, the behavior of the soliton inside the medium is 
governed by the combination of the field and matter char- 
acteristics. In the paper [15|, the combination T ~ Ct p 
was proposed, so that for longer pulses, one can obtain 
the effect at lower density. We study this relationship in 
Fig. E^b), where the probe pulse peak intensity is plotted 
versus the medium density (in units of Cq corresponding 
to the previously used Lorentz frequency wi = 10 11 s _1 ). 
One can see that at t p = 50 fs, the diode action covers a 
wide region of densities 0.7 < C/Co < 2.3 [B region, in 
terms of Fig. [9fa)]. At lower and higher densities, there is 
no directional asymmetry in pulse behavior (elastic and 
strongly inelastic conditions, respectively). If we take 
longer pulses (t p = 100 fs), region B shifts toward lower 
densities and becomes narrower, in accordance with the 
T combination. It is also important that the contrast 
between the curves in region B is not as impressive as 
for shorter pulses. However, we can say with confidence 
that the asymmetric propagation appears at lower densi- 
ties than previously due to the increase in pulse duration. 
In addition, a narrow C region appears at t p — 100 fs, 
while region A is located near the very vertical axis. In 
general, the T combination satisfactorily describes the 
effects of density and pulse duration on the location and 
width of the region of the diode action. 



CONCLUSION 



FIG. 9. (Color online) The dependence of the peak intensity 
of the probe pulse on (a) the duration of the pulses, and (b) 
the density of the medium (for two different durations). The 
value Co used for normalization corresponds to the Lorentz 
frequency lul = 10 11 s _1 . In figure (a) the value C = Co is 
adopted. 



controlled absorption and all-optical diode action. The 
results of calculations are depicted in Fig. M,&) . We mark 
out there several regions denoted with letters from A to 
D. In region A (very short pulses, t p < 30 fs) both FP 
and BP probe pulses appear at the output. Moreover, the 
total absorption of the first FP soliton is not observed 
in this case. Therefore, one can say that region A is 
the region of elastic collisions of the counter-propagating 
solitons, which is characteristic for short enough (high- 
intensity) pulses [Io|, OH. The region D (long pulses, 
t p > 150 fs) is the reverse case of strongly inelastic colli- 
sions, so that the soliton is annihilated independently of 
the propagation direction (l6| . Finally, the intermediate 
regions B and C are of interest for us. In region B (con- 
ditions considered above in this paper), the FP probe 
pulse is absorbed, while the BP one is transmitted. In 
region C the situation is reversed, so that one can change 
the direction of the diode action simply by changing the 



In conclusion, we studied the collisions of multiple 
counter-propagating SIT solitons within the framework 
of the semiclassical model of light interaction with a ho- 
mogeneously broadened two-level medium. It is shown 
that by means of inelastic asymmetric collision of two 
proper pulses, the medium can be prepared in such a 
state that it behaves as a diode transmitting the probe 
pulse or not, depending on propagation direction and 
peak intensity. This medium preparation is connected 
with the controlled light absorption in the medium which 
is limited by the process of saturation. By combining the 
characteristics of the medium and the pulse, one can find 
the proper conditions to observe the effects of asymmet- 
ric collisions. We should also notice that pulse retrieval 
was not achieved, so that the stored energy leaves the 
medium naturally as a consequence of fluorescence. 

One can expect that the effects of SIT solitons colli- 
sions can be observed in usual self-transparency media: 
solid dielectrics (including ruby where SIT was first dis- 
covered), vapors of alkali metals, molecular gases, and 
semiconductors. A review of the early experiments can 
be found in Ref. @. Obviously, the solid-state systems 
seem to be the most convenient for our aims. The con- 
ditions considered in this paper (very short and high- 
intensity pulses, short distances of propagation) are in 
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good agreement with the requirements of semiconduc- 
tor materials which possess short relaxation times and 
comparatively high transition dipole moments. The the- 
oretical and experimental results on self-induced trans- 
mission in semiconductors were reported, for example, in 
(l7l - l20l |. Another prospective material for SIT collision 
experiments is the collection of semiconductor quantum 
dots which can be considered as artificial two-level atoms 
with high dipole moments [U [HJ . 

The final remark deals with the problem of inhomo- 
geneous broadening when the resonant frequency of two- 
level atoms is distributed in a certain range. Recall that 
our calculations were performed for the case of a homoge- 



neously broadened medium, i.e., for broadening only due 
to the finite phenomenological relaxation Ti (equal for ev- 
ery atom). As a matter of fact, we are in the well-studied 
regime of SIT solitons which have essentially the same 
main properties in both cases of homogeneous and in- 
homogeneous broadening [s}. Moreover, the well-known 
area theorem (one of the attributes of SIT) is not valid in 
a strict sense for the homogeneously broadened medium 
, as it was confirmed recently by direct numerical sim- 
ulations [2j3|. Therefore, although the role of inhomoge- 
neous broadening is still to be studied, we believe that 
the main features of SIT pulse collisions should remain 
unchanged in that case. 
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